Elemental x-ray imaging using Zernike phase contrast
Shao Qi-Gang1, Chen Jian1, Wali Faiz1, Bao Yuan1, Wang Zhi-Li2, Zhu Pei-Ping3, †, , Tian Yang-Chao1, ‡, , Gao Kun1, ¶,
National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China
School of Electronic Science and Applied Physics, Hefei University of Technology, Hefei 230009, China
Institute of High Energy Physics, Chinese Academy of Science, Beijing 100049, China

 

† Corresponding author. E-mail: zhupp@ihep.ac.cn

‡ Corresponding author. E-mail: ychtian@ustc.edu.cn

¶ Corresponding author. E-mail: gaokun@ustc.edu.cn

Project supported by the National Basic Research Program of China (Grant No. 2012CB825801) and the National Natural Science Foundation of China (Grant Nos. 11505188, and 11305173).

Abstract
Abstract

We develop an element-specific x-ray microscopy method by using Zernike phase contrast imaging near absorption edges, where a real part of refractive index changes abruptly. In this method two phase contrast images are subtracted to obtain the target element: one is at the absorption edge of the target element and the other is near the absorption edge. The x-ray exposure required by this method is expected to be significantly lower than that of conventional absorption-based x-ray elemental imaging methods. Numerical calculations confirm the advantages of this highly efficient imaging method.

1. Introduction

X-rays have been employed as a quantitative mapping tool of elemental constituents since the last century.[1,2] At the absorption edge of an element, the complex x-ray atomic scattering factor experiences rapid and substantial change.[3,4] This elemental specific change has been exploited using x-ray absorption contrast imaging technique to obtain element specific image.[59] The technique typically requires two x-ray images: one is at x-ray energy slightly below the absorption edge energy and the other is at x-ray energy slightly above the absorption edge energy. After proper normalization of the two images in intensity, subtraction of the normalized images will yield an image of the target element. This elemental imaging technique is simple, straight forward, and fast, especially when a synchrotron x-ray source is used.[1013]

For high Z elements, however, the change of the real part of their atomic scattering factors across an absorption edge between 0.3 keV–2 keV can be larger than that of the corresponding imaginary part,[3,4] and we refer to the change as resonance peak. Figure 1 shows the components of the refractive index of four elements each as a function of x-ray energy.[4] Here, we underline that for each element, the change of the real part across an absorption edge is larger than that of the imaginary part. Effective use of the large change in the real part of the atomic scattering factor may therefore offer significant advantage in terms of elemental sensitivity, yield, and radiation dose.

Fig. 1. Plots of refractive-index decrement δ and imaginary part β of the refractive index versus energy for different metals (black for Ti, green for Mn, purple for Ni, and red for Zn). Curves are based on Ref. [4].

Here, we propose a new elemental imaging method by using Zernike phase contrast (the changes of both the real and the imaginary parts of the atomic scattering factor) instead of absorption contrast (solely the change of the imaginary part). In our proposed method, two phase-contrast images are collected by using Zernike phase contrast setup: one is at energy corresponding to the resonance peak and the other is close to the resonance peak. The two images are normalized and subtracted to obtain an image of the target element. Our proposed method can be readily performed on existing Zernike phase contrast x-ray microscopes that are available by using synchrotron and laboratory x-ray sources. In Zernike phase contrast x-ray microscopes, the sample is placed on the focus of the condenser, a zone plate is used as an objective, and a phase ring is placed in the back focal plane of the zone plate to achieve π/2 or 3π/2 phase shift.[1417] Only the working energy needs to be changed once in the system to perform the new elemental imaging experiment. Possible problems, such as different magnifications of the images at different energies, in the experiment could be solved just as in the absorption based method.[18]

2. Simulation and calculation

To demonstrate the effect of the proposed method, we simulate the imaging procedure of different metal particles within a model cell as illustrated in Fig. 2(a). The cell contains zinc (Zn) and gold (Au) nanoparticles each with a 20-nm diameter. The cell contains 80% water and 20% organic material with a composition of C22H10N2O5, and has an average density ρc = 1.43 g/cm3.[19] To test the element-sensitivity a protein cluster is placed near these particles. For the model protein structure we use the stoichiometric composition H48.6C32.9N8.9O8.9S0.3 with the density ρp = 1.35 g/cm3.[20] In the simulation, the setup is the same as that used in conventional Zernike phase contrast microscope.[1417] The resolution of the microscope is 20 nm, a gold-made phase ring is used, and the radius and width of phase ring are selected according to the condenser properties. Poisson noise is added to simulate the physical noise of x-ray photons. Since the phase shift of Zn changes dramatically near the resonance peak around 1022.3 eV,[3] we choose a suitable pair of energies, i.e., 1022.3 eV and 973.7 eV, as working energies. For comparison, we also calculate the required exposure of absorption methods. In the absorption method the working energies are 1022.3 eV and 1012.4 eV. Figures 2(b)2(d) show the simulation results of regular absorption transmission x-ray microscopy (TXM) images (Figs. 2(b) and 2(c)) and phase contrast image (Fig. 2(d)), in which Zn and Au particles are not clearly distinguished from each other. Figure 3 shows images after the subtraction by using the phase contrast data (Fig. 2(a)) and the absorption data (Figs. 2(b)2(f)), respectively. In subtraction images Zn nanoparticles stand out, while Au nanoparticles and the protein cluster disappear. In addition, we calculate the signal-to-noise ratios (SNRs) for the three Zn nanoparticles in Fig. 3 under different photon number densities (PNDs). The signal is an average value of 9× 9 pixels (the size of nanoparticles) and the noise is calculated by the standard deviation of the same size of background area. The data show that to achieve the same picture quality, the required exposure in phase contrast is only about one-third that required in the absorption method.

Fig. 2. (a) Profile of the sample, with the bulk protein and Zn and Au nanoparticles inside a cell; (b) absorption TXM image in the water window; (c) absorption TXM image at the energy E1 (1022.3 eV); (d) x-ray phase contrast image at the energy E1 (1022.3 eV). All images have PND = 4 × 1010 photons/μm2.
Fig. 3. Subtraction images: (a) phase contrast, and (b)–(f) absorption contrasts. The sample exposures correspond to PND: 4 × 1010, 4 × 1010, 8 × 1010, 1.2 × 1011, 1.6 × 1011, 2 × 1011 photons/μm2, respectively. The SNRs of the images are 12.6, 7.3, 9.9, 11.9, 13.8, and 15.4, respectively.

For in-depth analysis and quantifying the advantages of elemental imaging with using the phase contrast, we consider PND to achieve given SNR, as a criterion. The PND for conventional TXM images is calculated by[15]

where Is and Ib are the intensities of sample and background respectively, F is the area of the target feature, and ηzp is the first order efficiency of zone plate. Similarly, PND for the new elemental imaging method can be calculated as

where K = Ib2/Ib1 is the background normalization factor, and parameters with subscripts 1, 2 indicate parameters at corresponding working energies. For Zernike phase contrast system, the intensity downstream of the sample is given by[21]

and the intensity of the image of the matrix (background) is

where A0 is the complex amplitude on the incident plane, μs, μb, and μp are the absorption coefficients of sample, matrix, and phase ring, respectively, ηs, ηb, and ηp are the phase shift coefficients, t is the thickness of the sample, and tp the thickness of the phase ring.

3. Results

In order to show the performance improvement in terms of the dose, we calculate the PND ratio: nr = PND(absorption)/PND(phase). The model consists of a small Zn particle in a 5-μm cell. Working energies for the calculation are 1022.3 eV and 973.7 eV, 1022.3 eV and 992.9 eV, 1022.3 eV, and 1012.4 eV for the phase contrast and 1022.3 eV and 1012.4 eV for the absorption. The variations of PND ratio with nanoparticle size are shown in Fig. 4 correspondingly. Taking 20-nm Zn particle for example, the PND ratio of the absorption contrast imaging to phase contrast imaging can be up to 4.3, which means that dose can be reduced down to ∼ 1/4 that with the absorption method. We should also notice that the required exposure for the new method decreases as the difference between the two energies increases. This is because the difference in phase shift increases as one energy moves away from the absorption edge. However, this is not the case in the absorption method, thus we can further optimize the results by broadening the energy spacing in a certain range.

Fig. 4. Performance improvement in terms of exposure. The variations of ratio of the required x-ray exposure (PND) in absorption contrast to that in phase contrast with the size of the nanoparticles. Energies for absorption are 1022.3 eV and 1012.4 eV, and energies for phase contrast are indicated in the upper right corner of the figure.

We also perform similar calculations for other metals such as Cu, Fe, and Ti, data based on Ref. [4]. Results indicate that for small size particles (<100 nm), elemental imaging based on the phase contrast is more efficient, i.e., at constant SNR the PND for phase contrast is halved or even quartered compared with that in the absorption method.

In addition, we evaluate the performance of the new method in terms of the size of the particles, and the results show that our phase contrast technique is always more effective than the absorption method, as long as the object particle size is smaller than 100 nm. Our new method is suitable for the soft x-ray range. Also, we should mention that the data from Ref. [4] are not exactly accurate near absorption edges, but for metals they are close to the realistic ones. Since no molecular orbitals are involved, we will consider these situations in our following work.

4. Conclusions

Our proposed elemental x-ray phase contrast imaging method exploits the intrinsic phase changes of resonance peaks, which enables us to achieve low dose element-specific x-ray imaging. Despite several limitations, this new imaging method can be used in a wide range of biological applications, such as three-dimensional (3D) imaging of metal protein clusters, and even in material sciences, such as transition process of metal clusters in battery materials. Finally, our proposed method can serve as a good complement to XANES imaging and XFM techniques, and easily implement with current TXM microscopes.

Reference
1Engstrom A1946Acta Radiol. Suppl637
2Kennedy J MJacobsen CKirz JRarback HCinotti FThomlinson WRosser RSchidlovsky G 1985 J. Microsc-Oxford 138 321
3Henke B LGullikson E MDavis J C 1993 Atom. Data Nucl. Data 54 181
4http://henke.lbl.gov/optical_constants/
5Klems J H 1991 Phys. Rev. 43 2041
6Bernhard PMaul JBerg TWegelin FOtt USudek CSpiecker HMerchel SSchonhense G 2006 Phys. Rev. 74 075401
7Grew K NChu Y SYi JPeracchio A AIzzo J RHwu YCarlo F DChiu W K S 2010 J. Electrochem. Soc 157 B783
8Yin G CTang M TSong Y FChen F RLiang K SDuewer F WYun WKo C HShieh H D 2006 Appl. Phys. Lett 88 241115
9Fan J DJiang H D2012Acta Phys. Sin61218702(in Chinese)
10Sun LHuo YZhou CLiang J HZhang X ZXu Z JWang YWu Y Z 2015 Acta Phys. Sin 64 197502 (in Chinese)
11Chen JGao KWang Z LYun W BWu Z Y 2015 Chin. Phys. 24 086804
12Zhang KLi D EHong Y LZhu P PYuan Q XHuang W XGao KZhou H ZWu Z Y 2013 Chin. Phys. 22 076801
13Shahzad FSiddiqi S AIm M YAvallone AFischer PHussain ZSiddiqi IHellman F 2010 Chin. Phys. 19 037504
14Chen Y TChen T YYi JChu Y SLee W KWang C LKempson I MHwu YGajdosik VMargaritondo G 2011 Opt. Lett 36 1269
15Schneider G 1998 Ultramicroscopy 75 85
16Schmahl GRudolph DGuttmann PSchneider GThieme JNiemann B 1995 Rev. Sci. Instrum 66 1282
17Zhang KLi D EHong Y LZhu P PYuan Q XHuang W XGao KZhou H ZWu Z Y 2013 Chin. Phys. 22 076801
18Liu Y JMeirer FWang J YRequena GWilliams PNelson JMehta AAndrews J CPianetta P 2012 Anal. Bioanal. Chem 404 1297
19Wang Z LGao KChen JHong Y LGe XWang D JPan Z YZhu P PYun W BJacobsen CWu Z Y 2013 Biotechnol. Adv 31 387
20London R ARosen M DTrebes J E 1989 Appl. Opt 28 3397
21Duke P JMichette A G1990Modern Microscopies – Techniques and ApplicationsNew YorkSpringer59